New scalar field models and their defect solutions 
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In this work we introduce new scalar field models and study the defect solutions they may engen- 
der. The investigation is based on the deformation procedure, which greatly simplify the calcula- 
tions, leading us to new models together with the corresponding solutions in a simple and powerful 
manner. 
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Generalities. — Topological structures [l| are of gen- 
eral interest in physics because they may be used in a 
diversity of situations, ranging from high energy physics, 
where they can appear during phase transitions governed 
by spontaneous symmetry breaking [2j], going down to 
low energy physics, where they can contribute to pattern 
formation [3[ and to the construction of logical apparatus 
at the nanometric scale Q. 

In this work we focus on new models described by real 
scalar fields and their topological solutions. The inves- 
tigation is of current interest to a diversity of scenarios 
both in high energy physics and in condensed matter. In 
high energy physics, they may help understanding inter- 
esting problems related to topological defects; see, e.g., 
Refs. i"). 6]. In condensed matter, we can use the models 
to investigate interesting issues, for instance, to describe 
polarons in carbon nanotubes 0] and to study nonlin- 
ear excitations in Bose-Einstein condensates [H, Q . By 
addressing a much broader class of potentials we hope 
to provide analytical basis to the more and more exotic 
structures arising in the areas of physics previously men- 
tioned. 

Here we will take advantage of the deformation proce- 
dure introduced in Ref. 110] for real scalar fields in the 
standard case, dealing with systems defined in the real 
line. We consider models described by the Lagrange den- 
sities 



equation of motion 



c d = -d^d^~v{4>), 



(la) 

(ib) 



where x an d 4 1 are t w0 real scalar fields and U(x) an d 
V(4>) are given potentials, which specify each one of the 
two models. The approach requires that we introduce 
a function x = /(</>) j named deformation function, from 
which we link the model (|la|) with the deformed model 
(fib]) . This is done relating the two potentials U(x) and 
V((f>) in the very specific form 
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dx 2 



dU 

d x' 
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then the deformed model has finite energy static solution 
given by 4>{x) = f {x{ x ))i which obeys 



dx 2 



dV 



(4) 



The proof of the above result can be found in Ref. [T(| , 
so we omit it here. 

We notice that if the deformation function /(</>) 
has critical points, the deformation procedure works 
smoothly if and only if the potential U (f{4>)) has zeros of 
multiplicity at least two at those critical points. In this 
case, the critical points of /(</>) are the branching points 
of / (x), which is consequently a multivalued function; 
see, e.g., Refs. flll - [l3j for other issues concerning the de- 
formation procedure. 

New models. — We now proceed to one important 
point of the present work, which is to find new models 
and their defect solutions. To do this, we start with poly- 
nomial and irrational potentials, and we search for the 
corresponding kinks and/or lumps, which are the clas- 
sical static solutions of the corresponding equations of 
motion. We follow the standard methodology, taking as 
the starting model a family of models, described by the 
potential 



U(x) 



2\n 



(i-x 2 ) 
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Here we use dimensionless field and coordinates, and we 
take n as integer or half integer, such that n > 3/2. 

This is the new family of models which we consider, 
and for n — 2 we get to the standard fourth order power 
potential with spontaneous symmetry breaking. If we 
choose the center of the defect at the origin, we get to 
the kink (+) and antikink (— ) solutions ± Xn (x), where 
Xn{x) can be expressed in terms of the inverse of the 
Gauss Hypergeometric function, in a way such that 



The choice leads to the result that if the original model 
has finite energy static solution x( x ) which obeys the 
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In the case of n — 2 and n = 3, the above expression is 
easily invertible to give the explicit solutions 



and 



X2(x) = tanh(x) , 



X3(x) 



vi + x 2 



(7) 



(8) 



respectively. The case n — 2 is well known and n = 3 is 
nonstandard, but it was already considered in Ref. 

In general, it is very hard to write the explicit solution 
analytically, but the implicit plot can be given; see, e.g., 
Fig. 1 for the cases of n = 3/2, 5/2, and 3. To go further 
on this, we can write the potential in the form 



TT( \ 1 ( dWn 

U{X) = 2 {-W 



where 



W n ( X ) = 2F1 



n 3 
2' 2' 



X 2 I X- 
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Note that formula (9) is well defined on the reals for every 
X G (—00,00) only if n is an even integer. If n is integer 
but odd U is always real but W is purely imaginary for 
|x| > 1. Finally, if n is half-integer U and W in formula 
(9) are real only if x G [—1,1]. In any case, this leads to 
the energy of the static solution 



E n = \W n ( X = 1) - W n ( X = -1)| = 



(11) 

where -B(a, 6) stands for the Euler Beta function. For n 
even: n — 2m, with m 6 N, the corresponding superpo- 
tential pH|) . reduces to a polynomial in x 2 



W 2m ( X ) 



X 



2/V 



2k + 1 
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For n odd: n = 2m + 1, the expression for the superpo- 
tential (fTU|) can be expressed in the form 



II 



2m+l 



(x) 



(2m + 1)!! 
(2m + 2)!! 



arcsinx + X V 1 - X 2 P m (x 2 ) 

(13) 

where P m (x 2 ) 1S a polynomial of degree m in % 2 , thus of 
degree 2m in the field x- 

Deformed models. — In the previous section, we 
have constructed new family of models. Thus, let us now 
focus on the main point of this work, which concerns the 
construction of new models, with the use of the deforma- 
tion procedure. We use the deformation function [ll| 



f{4>) = cos(a arccos cj> — mm), 



(14) 



where a is a real constant, and m can be integer or half 
integer. We deform the potential (|S|), which depends on 
n, with the deformation function (|14[) . which depends on 



// 


: \ 


f -1.0 -0.5 




0.5 1.0 ^\ 






1.0 






0.5 






-2 -1 / 
-0.5 




1 2 



FIG. 1. In the upper panel: plots of the potential for n = 3/2 
(large dashing), 5/2 (small dashing), and 3 (continuous). In 
the lower panel: plots of the respective kinklike solutions. 



a and m. The parameter m leads to two distinct families 
of models: for m integer, the deformed potential can be 
written in the form 



1 



(l-^)sin 2 "- 2 ( 



a arccos 



4>) (15) 



and, for m half integer we get 



V 

1 r 



11 .a 

cos 



W>) = 2^2(1 -0 2 )cos 



2n— 2 



(a arccos 0) (16) 



The new models (|15|) and (|16p do not depend on m any- 
more, since we have already used the fact that m is in- 
teger or half integer. Moreover, we will see below that 
using integer or half integer values for the parameter a 
leads the number of minima of the new potentials to be 
fixed at will. 

The above models identify families of potentials which 
present static solutions which can be constructed with 
the help of the general expression 



i(x) + m 7r 



where 



^(x) = arccos(xn(^)), 



(17) 



(18) 



with Xn {%) being defect solution of the model discussed 
in the previous section, and 9 n (x) € [0, tt]. 

Since we have a diversity of models and solutions, let 
us now illustrate the investigation focusing our attention 
on some specific cases. 

New models for a integer.— Here we investigate 
families of models for the specific cases of a and n being 
integers. We first consider the case of V^°, as written in 
flT5"l) . We have two possibilities: 



3 



• for a odd: 



2a 

• for a even: 



(a+l)/2 / 2 V 2n-2 

b^r n 



(19) 



_-^- 2 a-^rjg(i-|.J 

(20) 

where Z° = cos (^pTr)- The sine potential can also be 
written in terms of the Chebyshev polynomials of second 
kind, in the form 



2n-2 



where 



sin((a + 1) arccos< 
sin(arccos 9) 



(21) 



(22) 



The explicit form of the potentials are given by, for a — 
1,2,3 and n = 2 



(23) 
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FIG. 2. Plots of V s "; a for some n and a. 
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and for 



K^W=2 4 (1-^) 
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(24) 



(25) 



(26) 



(27) 
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We notice that for n even, there are two classes of models: 
for a odd they are (j) 4 — like potentials (no zero at the 
origin) and for a even they are (f> e — like models (having 
a zero at the origin). And for n odd, there are also two 
classes of models: for a odd they are inverted 4> 4 — like 
potentials (no zero at the origin) and for a even they are 
inverted </> 6 — like models (having a zero at the origin). 
To illustrate the results, in Fig. 2 we plot some of the 
potentials as a function of the scalar field cf>. 

The static solutions of the potentials V™^ a (4>) are given 

by 



4>k(x) = cos 



9 n { X ) , (fc-l)7T 



(29) 



where fc is integer, which produces different solutions only 
if 1 < k < 2a; recall that (9 n (x) is given by (JTSJ. We 
notice that, for n even, the potentials of the sine family 
are non negative, and all the zeros are minima. Then, 
all the static solutions given by Eq. (|2T)| are topological, 
since they interpolate consecutive minima. In the case for 
n odd, the potentials of the sine family are non negative 
only for \(f>\ < 1, and the zeros at <j> = ±1 are not minima. 
Then, all the static solutions given by Eq. (|29| which 
interpolate consecutive local minima are topological, and 
the others, for k = a, 2a, which involves the zeros at <j> = 
±1 and the neighbor local minimum, are non topological. 

Let us now consider the family of models V^.™i, a given 
by (|T6|). for the specific cases of a and n being integers. 
Here we also have two possibilities: 



• for a odd: 



Kos"(0) 



• for a even: 



(a-l)/2 

2n - 2 {l-4> 2 ) J] 
J'=l 



(30) 



(31) 



cos (^ji7r). The cosine potentials V™ f{4>) 
can also be written in terms of Chebyshev polynomials 



where 



4 



of the first kind 



where 



T a {9) = cos(aarccos( 



(32) 



(33) 



The explicit form of the potentials are given by, for 
a = 1, 2, 3 and n = 2 





(34) 



(35) 




and for n = 3 



(36) 



(37) 




K 3 o f(0)^(W 2 ) Q-^ 2 



128 



4 (1-0 2 ) 



(38) 



We notice that there are two classes of models: for a 
even they are inverted 4 — like potentials (no zero at the 
origin) and for a odd they are inverted <fP— like models 
(having a zero at the origin). To illustrate the results, in 
Fig. 3 we plot some potentials as a function of the scalar 
field cj). 

The static solutions of the V££ (</>), are given by 



<J>k(x) 



t (x) (2fc-l)7T 



2 a 



(40) 



where k is integer, which produces different solutions only 
if 1 < k < 2a; recall that 9 n (x) is given by fl8]). We 
notice that, the potentials of the cosine family are non 
negative only for — 1 < (f> < 1, and the zeros at (f> — 
±1 are not minima. Then, for n integer, all the static 
solutions given by Eq. (|4"0"|) which interpolate consecutive 
minima are topological, and the others, for k — a, 2a, 
which involves the zeros at <fi = ±1 and the neighbor 
minima, are non topological. 

New models for a half integer. — In this case, the 
two families of models are field reflection of one another, 
that is, V^f((j)) = K,"; a (— </>). Thus, they have no impor- 
tant difference, and we will examine explicitly only one of 



them. The potentials of the sine family can be expressed 
in the form 



(39) C(« 



1 



2n-2 



2 n a< 



(1-0 2 )(1- 



(41) 

where Zj = cos ( \ a tt); and n is an integer. They can 
also be given in terms of Chebyshev polynomials of the 
first kind 



1 



(42) 



where T a (9) is given by Eq. (|3"3")l . The explicit form of the 
potentials are given by, for a = 1/2, 3/2, 5/2 and n = 2 



^ /2 («/>) = (l-^) (1- 



K?a /2 W= g(W 2 ) (I"*) ( 1 + 2 ^) 2 



(43) 
(44) 



Kl 5/2 W = ^ (1 - 2 )(1 - 0)(1 - 20 - 40 2 ) 2 , (45) 



and for n = 3 



K-; 1/2 W) = ^(i-0 2 ) a 



Kin 2 (0) = I g(l-0 2 )(l-0) 2 (1+20) 4 



(46) 
(47) 
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V^'\<t>) = ^ (1 - 4> 2 ){l ~ ^) 2 (l - 20 - 40 2 ) 4 , (48) 

In Fig. 4, we plot some of the potentials as a function 
of the scalar field (f>. Also, the static solutions in this case 
are given by 



4>k(x) = cos 



l n {x) , (fc-l)TT 



where fc is integer, which produces different solutions only 
if 1 < k < 2a; recall that 6 n (x) is given by (fT5)) . 

We notice that these potentials are non negative only 
for <f> > — 1, and the zero at 4> = —1 is not a mini- 
mum. Thus, all the static solutions given by Eq. (|49)) 
which interpolate between consecutive minima are topo- 
logical, and the other, for k = a + 1/2, which involves 
the zero at <j> = —1 and the neighbor minimum, is non- 
topological. 



not non negative. Nevertheless, we can follow the lines 
of [lj] and introduce superpotential, for both topological 
and for non topological solutions. In the case of the sine 
families, for a integer or semi-integer, the superpotentials 
can be written in terms of Chebyshev polynomials if we 
fix the integer n, for instance: 



(49) W 2,a u , = (a 2 (l - </> 2 ) - 2)T a {<P) - 2a0(l - 4> 2 )Ug-i{4>) 



a 2 (a 2 -4) 



for a ^ 2, and 



(51a) 



_ <jVj - 4> 2 + arcsin^ y/l - 4> 2 
4a 8a 2 (a 2 - 1) 



x [2a0r 2a (0) - (1 - 2a 2 (l - 2 ))f/ 2a _ 1 (c^)] 

(51b) 

Also, if we fix a we can get other expressions for n integer, 
for instance: 
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FIG. 4. Plots of V™^ a for some n and a. 



Superpotential and energy. In general, when 
the potential is non negative, it is possible to introduce 
superpotential W = W(<j>) such that 



v{4>) 



1 fdW 

2 V"#" 



(50) 



This is the case for the sine family of potentials with a in- 
teger. However, In the others cases, for the cosine family 
of potential with a integer, and for the sine and cosine 
families of potential with a half integer, the potentials are 



(52b) 

where F\ is the Appell function of first type. Also, in the 
case of the cosine families, for a integer or semi-integer, 
the expressions for the superpotentials are similar to the 
corresponding expressions for W a \ n , for instance: 



i 



a 2 (a 2 - 4) 



[(a 2 (l — <p 2 ) — 2) sin(a arccos<? 



+2a(f>\/ 1 — 4> 2 cos(a arccos</>)] 
for a ^ 2, and 



(53a) 



n+l /l 1 3 

2— ^ 1 U ; -2' 1 - n; 2 ; ^ 2 ^ 

(53b) 



( i yi+iQn-3 j.n 



(n + 2) 



, n 1 n 

" 1 2 + 5 ~2' 1_n; 2 +2 



77- I 77- 

-3(n + 2)fi ( 2 ; ~2'~ n; 2 +1 



(53c) 



In general, the superpotential simplify the calculations. 
Particularly, for the energy associated with the corre- 
sponding static solutions we can write, for the kinklike 
solutions, 

E n ^ k {<jy k ) = \W n ' a (Z k ) - W n ' a (Z k+1 )\ , (54) 
and for the lumplike solutions [lij 

E n ^\cp k ) = 2 \W n > a {Z k ) - W n ' a {4>{x = 0))| . (55) 
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The general superpotential, for arbitrary n and a, is 
very hard to be calculated. However, in order to calculate 
the energy, we can circumvent this problem following the 
alternative route: the kinklike solution Xn( x )i Eq.©, can 
be approximated by 



Xn{x) = ±tanh(a„x) 



(56) 



where the width of the solution is controlled by the in- 
verse of 



3tt 1 



(57) 



This is an approximation, but it leads to an error less 
than two percent for the energy of the corresponding so- 
lution. Also, we can write the energy as 



En 



dx 



dx 



(58) 



With this last expression, we use ([56]) to get that E n = 
4a„/3. Now, comparing the approximated energy value 
with those from dill) , we obtain a good fit by means 
of (|57|) . Particularly, for n = 3/2,5/2 and 3 we obtain 
a 3/2 = 1.052, a 5/2 = 0.912, and a 3 = 0.866. For the 



models described by the potential ([5]), the energy of the 
defect solution decreases with increasing n, but the corre- 
sponding width increases. Finally, for n^2, the energies 
of the static solutions given by Eqs. (f2"9"|). (HOT) , and (|4"9"j) . 
are obtained straightforwardly from the energy for n = 2, 
with the answer E n > a > k = a n E 2 - a - k . 

Final comments. — In this work, we have introduced 
and solved two distinct new families of models. The first 
family is defined in ([5]), and it is controlled by the single 
parameter n, integer or half integer. The second family 
is introduced by the deformation function (|14|) , which 
depends of two new parameters, a and to. With this 
deformation we could find new models, depending on m 
being integer or half integer, and on a and n. In this 
case, we have obtained a diversity of models and their 
respective defect solutions. We have also found a way of 
presenting the energy of the defect solutions, both exactly 
or approximately, when the exact result is hard to be 
obtained. As one knows, we can also investigate linear 
stability of the scalar field solutions. Although this issue 
is out of the scope of the present work, it can be done 
following the lines of [Til ]. 
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